Search for: All records

Abstract The nonlinear mechanical responses of rocks and soils to seismic waves play an important role in earthquake physics, influencing ground motion from source to site. Continuous geophysical monitoring, such as ambient noise interferometry, has revealed co‐seismic wave speed reductions extending tens of kilometers from earthquake sources. However, the mechanisms governing these changes remain challenging to model, especially at regional scales. Using a nonlinear damage model constrained by laboratory experiments, we develop and apply an open‐source 3D discontinuous Galerkin method to simulate regional co‐seismic wave speed changes during the 2015 M_w7.8 Gorkha earthquake. We find pronounced spatial variations of co‐seismic wave speed reduction, ranging from <0.01% to >50%, particularly close to the source and within the Kathmandu Basin, while disagreement with observations remains. The most significant reduction occurs within the sedimentary basin and varies with basin depths, whereas wave speed reductions correlate with the fault slip distribution near the source. By comparing ground motions from simulations with elastic, viscoelastic, elastoplastic, and nonlinear damage rheologies, we demonstrate that the nonlinear damage model effectively captures low‐frequency ground motion amplification due to strain‐dependent wave speed reductions in soft sediments. We verify the accuracy of our approach through comparisons with analytical solutions and assess its scalability on high‐performance computing systems. The model shows near‐linear strong and weak scaling up to 2,048 nodes, enabling efficient large‐scale simulations. Our findings provide a physics‐based framework to quantify nonlinear earthquake effects and emphasize the importance of damage‐induced wave speed variations for seismic hazard assessment and ground motion predictions. 

SUMMARY The yield surfaces of rocks keep evolving beyond the initial yield stress owing to the damage accumulation and porosity change during brittle deformation. Using a poroelastic damage rheology model, we demonstrate that the measure of coupling between the yield surface change and accumulated damage is correlated with strain localization and the Kaiser effect. Constant or minor yield surface change is associated with strong strain localization, as seen in low-porosity crystalline rocks. In contrast, strong coupling between damage growth and the yield surface leads to distributed deformation, as seen in high-porosity rocks. Assuming that during brittle deformation damage occurs primarily in the form of microcracks, we propose that the measured acoustic emission (AE) in rock samples correlates with the damage accumulation. This allows quantifying the Kaiser effect under cyclic loading by matching between the onset of AE and the onset of damage growth. The ratio of the stress at the onset of AE to the peak stress of the previous loading cycle, or Felicity Ratio (FR), is calculated for different model parameters. The results of the simulation show that FR gradually decreases in the case of weak coupling between yield surface and damage growth. For a strong damage-related coupling promoting significant yield surface change, the FR remains close to one and decreases only towards the failure. The model predicts that a steep decrease in FR is associated with a transition between distributed and localized modes of failure. By linking the evolving yield surface to strain localization patterns and the Kaiser effect, the poroelastic damage rheology model provides a new quantitative tool to study failure modes of brittle rocks. 

Abstract We employ direct numerical simulations in order to analyze the role of double‐diffusive salt fingering in halite precipitation from hypersaline lakes. Guided by field observations from the Dead Sea, which represents the only modern deep stratified lake that precipitates halite under hydrological crisis, we consider a saturated layer of warm, salty brine (epilimnion) overlying a layer of colder, less salty brine (hypolimnion) that is also saturated. The double‐diffusive instability originating in the metalimnion gives rise to an asymmetrical pattern of less salty ascending fingers, accompanied by descending salt fingers that lose heat as they propagate through the metalimnion. The net result is a strong, downward salinity flux that leads to the undersaturation of the epilimnion, while the hypolimnion becomes oversaturated and precipitates halite. These interfacial processes within deep, hypersaline water columns in warm and dry regions suggest a potential route toward the formation of thick halite layers found in the geological record.